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PARTICLE ACCELERATION IN THE COLLIDING WINDS BINARY WR140 

J. M. Pittard, 1 and S. M. Dougherty 2 
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RESUMEN 

Favor de proporcionar un resumen en espanol. If you cannot provide a Spanish abstract, the 
editors will do this. Massive WR+O star systems produce high-temperature, shock-heated plasma where 
the wind of the WR star and that of its binary companion collide - the wind-collision region (WCR). The 
WCR is a source of thermal (e.g. hard X-rays) and non-thermal (e.g. synchrotron) emission, the latter arising 
from electrons and ions accelerated to relativistic energies. These colliding wind binaries (CWBs) provide an 
excellent laboratory for the study of particle acceleration at higher mass, photon and magnetic energy densities 
than exist in SNRs. Recent models of the non-thermal (NT) emission from WR 140 have provided insight into 
this process. 

ABSTRACT 

Massive WR+O star systems produce high-temperature, shock-heated plasma where the wind of the WR star 
and that of its binary companion collide - the wind-collision region (WCR) . The WCR is a source of thermal 
(e.g. hard X-rays) and non-thermal (e.g. synchrotron) emission, the latter arising from electrons and ions 
accelerated to relativistic energies. These colliding wind binaries (CWBs) provide an excellent laboratory for 
the study of particle acceleration at higher mass, photon and magnetic energy densities than exist in SNRs. 
Recent models of the non-thermal (NT) emission from WR 140 have provided insight into this process. 
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1. INTRODUCTION 

WR 140 (HD 193793) is the archetype of long- 
period CWB systems. It consists of a WC7 
star and an 04-5 star in a highly eccentric orbit 
(e w 0.88), and exhibits dramatic variations in 
its emission from n e ar-IR to radio wavel engths 
((Williams et all Il990t IWhite fc Becken Il995l ) and 
also at X-ray energies dZhekov fc Skinner! l200tj 



provided a more accurate representation of the 



thermal and NT e mission ((Dougherty et al. 120031 : 



Pollock. Corcoran fc Stevens! 120021: IPollock et all 
2005! ) during its 7.9-year orbit. The variability 
appears to be linked to the WCR, which experiences 
significant changes as the stellar separation varies 
between 2 and 30 AU. The orbit modulation of the 
synchrotron flux has yet to be understood, but is 
likely due to a number of mechanisms. The changing 
free-free opacity along the line of sight through the 
extended stellar winds is certain to play a role, as 
will the strong inverse Compton cooling of the NT 
electrons. In addition, the intrinisic synchrotron 
luminosity and the spectral index of the NT electron 
energy distribution may alter. Recently developed 
models of CWBs that are based on hydrodynamical 
simulations of the stellar winds and the WCR have 
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Pittard et al. 20_06 ,; ), and are a first step towards 
distentangling these mcchanims and ultimately 
understanding the acceleration processes in detail. 

Due to the non-unique solutions which can arise 
from these models, it is essential that observations 
across a broad energy spectrum are available to pro- 
vide as many contraints as possible. Fortunately, 
recent observations of WR 140 with the VLBA en- 
able a full orbit definition, including inclination, 
along with robust dista nce and luminosity estimates 
((Dougherty et al.ll2005r ). Together, these parameters 
provide essential constraints that are currently un- 
available for any other wide CWB. However, thermal 
X-ray and NT X-ray and 7-ray observations are also 
critical to establish some of the model parameters. 

2. MODELLING THE NON-THERMAL 
EMISSION 

Wc have applied our newly developed radiative 
transfer models of CWBs to WR140 in order to 
investigate the emission and absorption processes 
which act to go vern the radio variations. Fu ll details 
can be found in Pittard fc Dougherty ( 2006! ) . Rather 
than attempting to model the radio lightcurve in one 
step, our first aim was to obtain good fits to the ra- 
dio data at phase 0.837. This phase was chosen on 
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Fig. 1. Model fits to radio data of WR140 at orbital 
phase 0.837. In model A 77 = 0.22, while model B has 
77 = 0.02. The latter is preferred (see text for details). 

the basis that the observed radio emission is close 
to maximum, that an X-ray spectrum exists, and 
that orbital-induced curvature of the WCR is negligi- 
ble. The thermal X-ray flux is used to obtain a fam- 
ily of solutions with varying WR and O-star mass- 
loss rates as a function of wind momentum ratio, 
77. With these constraints, fits to the radio emission 
(see Fig.[I| lead us to favour models with a relatively 
small value of 77 (e.g. 0.02), on the basis that the re- 
quired line-of-sigh t angle is more consiste nt with the 
orbital solution of lDoughertv et ah (2005). Unfortu- 
nately, the VLBA images do not have the sensitivity 
to directly constrain the value of 77 due to the rapid 
decline in the surface brightness of the WCR with 
off-axis distance. In these models, free-free absorp- 
tion is responsible for the turnover between 1.6 and 
5 GHz - the Razin effect can be ruled out because it 
places an unacceptably large fraction of energy into 
NT electrons. The post-shock B-field at the apex of 
the WCR is estimated to be about 1 G, while some- 
what less than 1% of the total available energy is 
transferred to NT electrons. 

A key finding is that the slope of the energy spec- 
trum of the NT electrons, p, is flatter than the canon- 
ical value anticipated for diffusive shock acceleration 
(i.e. p < 2). Several mechanisms can produce such 
distributions. The most likely are re-a cceleration 
at m ultiple shocks within the WCR (e.g. ISchneider 
19931 ). and second-order Fermi acceleration result- 
ing from magnetic scattering off turbulent cells (e 
Scott fc Chevalierlll975l : iDolginov fc Silant'evlll99 
Both of these processes could occur together if the 
WCR is highly structured, as occurs naturally when 
clumps within the winds impact the WCR and in- 
troduce vorticity (jPittardl 120071 ) . Magnetic recon- 



nection, perhaps through resistive MHD, is another 
possibility. Reconnection probably occurs through- 
out the volume of the turbulent WCR, and not just 
at a hypothetical contact discontinuity, and may pro- 
vide additional energy for generating and maintain- 
ing the magnetic fluctuations which drive stochastic 
acceleration. 

A remaining problem is that the B-ficld and p 
are somewhat ill constrained by the radio data alone. 
However, their degeneracy is broken at 7-ray ener- 
gies, and thus tighter constraints can be made if fu- 
ture observations with GLAST are utilized. While 
we generally predict lower 7-ray fluxes than previous 
works, detectability should not be an issue. However, 
we conclude that WR 140 is unlikely to be detected 
as a TcV source with VERITAS-4, though it may 
be brighter at phases closer to periastron. Since the 
high stellar photon fluxes prevent the acceleration of 
electrons beyond Lorcntz factors 7 > 10 5 — 10 6 , TeV 
emission from CWB systems will provide unambigu- 
ous evidence of pion-decay emission from accelerated 
ions. 

3. FUTURE DIRECTIONS 
While these models have provided key insights 
into the particle acceleration process (es) occuring in 
WR140, much work remains. In forthcoming work 
we will apply our model to other orbital phases, and 
will investigate the effects of particle acceleration 
within the WCR. 
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